Compaction dynamics in ductile granular media 
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Ductile compaction is common in many natural systems, but the temporal evolution of such 
systems is rarely studied. We observe surprising oscillations in the weight measured at the bottom 
of a self-compacting ensemble of ductile grains. The oscillations develop during the first ten hours 
of the experiment, and usually persist through the length of an experiment (one week). The weight 
oscillations are connected to the grain-wall contacts, and are directly correlated with the observed 
strain evolution and the dynamics of grain-wall contacts during the compaction. Here, we present 
the experimental results and characteristic time constants of the system, and discuss possible reasons 
for the measured weight oscillations. 
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I. INTRODUCTION 

The stress distribution in dry granular media have been 
studied for more than a century. The German engineer 
Janssen studied the apparent weight at the bottom of a 
silo as function of its filling height Janssen found 
that the pressure at the bottom of a container of granu- 
lar material increases linearly with small filling heights, 
but approaches a constant level exponentially slowly for 
large filling heights. That the measured weight at the 
bottom is less than the total weight of grains is referred 
to as a screening effect. It is well known that the screen- 
ing effect is due to the grain-wall friction and how the 
stress distributes in a granular ensemble Q- Janssen's 
mathematical expression for this, the Janssen law, com- 
pares surprisingly well to experiments 0, , in spite of 
its crude assumptions regarding friction and stress dis- 
tribution p|. Over the last decade, various aspects of 
the stress distribution in static granular media have been 
studied. Experiments have shown that the stress distri- 
bution is sensitive to the deposition history |(| , the shape 
and size distribution of grains Q , elastic properties of the 
base [8| and grains |9|, and that an exponential size dis- 
tribution of forces is found at the bottom of a container 
for forces larger than the average 0, [n| . 

The stress distribution in dynamic systems has been 
investigated in pushed columns of granular media inside 
a cylinder 0, 0, ll3L If"H by measuring the normal force 
at the bottom for constant driving velocities. At small 
velocities, the measured force has a stick-slip behavior 
[l4l HE} that is related to aging of the grain-wall friction 
due to capillary condensation and shear strengthening 
of the contacts at the walls . These dynamic systems 
consist of elastic particles, and the time dependence stud- 
ied relate to other properties than the particle rheology. 
In Nature, and in many technological processes, slowly 
compacted or sheared systems are dominated by the de- 
formation of particles. The time dependence in these 
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systems is mainly given by the plastic properties of the 
grains. 

Here, the results from experiments on granular media 
consisting of plastically deforming grains in a cylinder 
are presented. This system deformed slowly under its 
own weight, compacting 10% in a week, while the nor- 
mal force at the bottom (the apparent mass |16|) was 
measured. The initial expectation was that the system 
would show a granular Janssen type stress distribution 
in the initial stage, but that due to the viscous rheol- 
ogy of the grains a stress distribution close to the hy- 
drostatic would develop. Thus, the apparent mass was 
expected to increase. Instead, the apparent mass devel- 
oped unexpected (non-harmonic) oscillations, resembling 
the stick-slip behavior observed in hard granular media 
pH Il5j , except that it decreased initially, and the time 
scale of a "slip" could be up to an hour. No overall in- 
crease was observed in the apparent mass after a week 
of compaction. The strain development and wall contact 
dynamics were also studied during the compaction, both 
showing behavior related to the weight oscillations. The 
wall interaction between grains and cylinder was varied 
significantly in a few experiments, and proved crucial to 
the oscillations, as these disappeared when wall friction 
was increased or decreased. The experiments and results 
are described in the following two sections, while some 
proposed mechanisms for the oscillations are discussed 
in section HVI 

We propose that the observed oscillations are due to 
two competing effects: Grain-wall interaction opposing 
motion and the slow flow and relaxation of the grains 
inducing motion. 



II. EXPERIMENTS 

We performed 30 experiments in which an ensemble of 
N deformable grains were left to compact in a Plcxiglas 
cylinder of diameter D. The system was studied in sev- 
eral ways, but mainly by measuring the apparent mass 
m a at the bottom of the cylinder in order to follow the 
overall evolution of the stress distribution in the com- 
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FIG. 1: (a) Schematic illustration of the setup. Ductile grains 
were filled to a height h in a cylinder and left to compact while 
the apparent mass m a at the bottom was measured, (b) A 
typical recording of the apparent mass (shown normalized by 
the total mass m of the grains) as a function of time. 



pacting ensemble. A Mettler PM4800 balance was used 
to measure the apparent mass. This balance operates by 
an induction mechanism that keeps the vertical position 
of the measurement area constant 0]> and thus does 
not affect the compaction procedure. The weight was 
measured to a precision of 0.03 g, and was typically a 
fraction (2-3)T0~ 4 of the total mass of the grains. The 
cylindrical container was mounted outside of the mea- 
surement area. Spherical grains were prepared manually 
from Play-Doh (Hasbro International Inc., UK) to a di- 
ameter d = (8.8 ± 0.2) mm, and poured into the cylinder 
approximately ten at a time. The initial packing fractions 
were in the range c = 0.5-0.6. The material is viscous 
[l8| over a large range of strain rates, e = (10 _2 -10~ 6 ) 
s _1 , with a viscosity of n = 3 ■ 10 5 Pas. A schematic 
illustration of the setup is shown in Fig. QJa) along with 
the typical result of the observed weight as a function 
of time. The measured apparent mass m a presented in 
Fig. mb) has been normalized by the total mass m of 
the grains in the cylinder. The apparent mass was found 
to oscillate in a quasi periodic manner. The period de- 
pended on details of the packing, and could increase or 
decrease slightly over the duration of each experiment. 
The filling height h(0) at t = was varied between 1-4 
times the cylinder diameter, and the cylinder diameter 
was varied between 3.4 and 15 times the grain diameter. 

In two experiments the total height h(t) of compact- 
ing (granular) ensembles were measured using two dif- 
ferent setups: A camera was used in one experiment to 
take pictures of the compaction process at various times. 
Image analysis was then used to extract the height of 
the ensemble based on the position of the uppermost 
6 grains, to a resolution of 46 /im. In another experi- 
ment, the height in the middle of the granular column 
was recorded by the use of a laser and two mirrors. A 
small, light weight piston was placed on top of the cen- 
tral grains and allowed to move only vertically. A small 
mirror was hinged onto the top of the piston, its lower 
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FIG. 2: Illustration of the experimental setup for strain mea- 
surement by the use of mirrors and laser; A balance (a) 
recorded the apparent mass m a at the bottom of the cylin- 
der. The height of the packing was measured as a function 
of time by a laser (b) beam that was reflected in a small and 
a large mirror (c), onto a point on the floor x. The position 
x moved to the left as the angle <f> between the small mirror 
and the horizontal plane was reduced, following the reduction 
of the height h of the compacting grains. The piston rested 
on grains that did not touch the walls, thus the strain was 
measured along a central axis. 



end resting on the upper cylinder rim. As the grains 
compacted, the mirror was tilted, reducing its angle <p 
to the horizontal plane. A laser beam was reflected in 
the small mirror, and again in another, larger, mirror so 
that the beam was visible as a point on the floor. The 
position of this point was recorded manually with time, 
and the height of the granular ensemble calculated to a 
precision of 3 /im. The piston was positioned along the 
central axis of the container, and followed the motion of 
the internal grains that did not touch the wall. Figure [2| 
illustrates the second strain measurement method. From 
the measurements of the total height the global strain, e, 
was found as e = 1 — h(t)/h(0). 

The dynamics of the grain contacts at the cylinder wall 
was studied using a camera (AstroCam, Capella, LSR 
Life Science Resources, UK) in one experiment. The 
camera had a spatial resolution of 2000 x 3000 square pix- 
els, and 14 bit intensity resolution. The contrast between 
the intensity inside and outside of a contact area was 
within an 8 bit subset of the 14 bit dynamic range. The 
rim of a contact was established within two pixels with 
the spatial and intensity resolutions as specified. The un- 
certainty that one extra rim of pixels introduced to the 
area of a contact could be as high as 20% for the smallest 
contact areas. The precision of the center of mass posi- 
tion was, however, much better, as it does not depend on 
the exact choice of thresholding for the contact area. 

The cylinder containing the ductile ensemble was 
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placed in front of two mirrors which were set in an angle 
of 72° to each other. The cylinder was reflected twice 
in each mirror, thus the camera view held five versions 
of the cylinder (I— V), capturing it from all sides. The 
grains' contacts to the wall were literally highlighted by 
shining light down the Plexiglas wall of the cylinder. The 
light only reflected out of the wall in areas where the dif- 
ference in refraction indices were smaller than that be- 
tween Plexiglas and air, thus the contacts between grains 
and wall were bright in contrast to the pore space. Fig- 
ure EUa) illustrates the setup. Each of the five cylinder 
images I-V (see Fig. |3J) was then 'unwrapped' and 
scaled according to the geometry of the setup, then put 
together to form a continuous image of the surface area of 
the cylinder. An example of the resultant image is shown 
in Fig. Elb). The spatial resolution in these images were 
160 /im. Images were recorded every 10 or 20 minutes 
for two days in order to capture several oscillations. 

A total of 90 contacts were recovered, and 79 of these 
were used in the analysis. The remaining 11 contacts 
were discarded because of some mismatch of their area 
across boundaries between cylinder images. 

An increase of the contact area of 70% was observed 
during the 60 hours that images were recorded, 60% dur- 
ing the first 20 hours of compaction, and 10% in the time 
interval t S [20, 60] h. A contact diameter was defined as 
for each contact area A, and found as a function 
of time. The average contact diameter, d c , was found by 
first taking the average value of each contact diameter 
over the series of time steps in t £ [20, 60], and then find 
the average of this set, d c — 2.66 ± 0.02 mm. 



III. RESULTS 

The typical behavior of the apparent mass m a in an 
experiment is as follows: At time t = all grains have 
been poured into the cylinder. The apparent mass in- 
creases slightly over a period of a few minutes, reaches 
its maximum (often a global maximum) and then starts 
to decrease. Weight oscillations mostly initiate during 
this decrease. When oscillations have developed, their 
minima decrease toward a global minimum of m a , be- 
fore they increase slowly toward a plateau. The plateau 
varies between experiments in the range 45%-88% of the 
total mass m of the grains, but are mostly in the range 
60%-80%. 

Figure 0] illustrates the definition of the periods, in- 
tervals, and amplitude of an oscillation, which will be 
referred to in the following. The period At of one oscil- 
lation is defined as the time between peaks i and (i + 1). 
This period can be further divided into intervals td and 
£j of overall decrease and increase, respectively, of the 
apparent mass. The point of minimum apparent mass 
between peaks % and (i + 1) marks the transition between 
the regions td and ij, see Fig. Etb). The amplitude Aa 
of one oscillation is the change in normalized apparent 
mass m a /m during ti. 




FIG. 3: (a) Schematic drawing of the setup for the measure- 
ment of contact areas at the wall of the cylinder. Two mir- 
rors in an angle 72° to each other reflect the cylinder surface 
and the total area can be extracted. Light emitting diodes 
were fitted into the top of the cylinder wall to enhance the 
contrast between contact regions (white) and regions of no 
contact (gray), (b) The unwrapped |Tsf| surface after image 
treatment. Each of the five (I-V) cylinder images is scaled 
and unwrapped before they are fitted in overlapping regions. 
The match is only at the cylinder surface, which is why the 
internal regions seems mismatched in some places. 



The weight oscillations initially have small amplitudes, 
Aa, which increase toward a maximum after typically 3- 
16 oscillations. The amplitudes reduce somewhat after 
this maximum value; In some experiments they nearly 
disappear after 100 hours, while in others they are still 
at their maximum value after 200 hours. The period At 
of an oscillation also tends to increase initially, and then 
stabilize at a constant value after typically 17-80 hours. 
In a few cases the period only stabilized after 150 hours, 
or not at all in the time span of the particular experiment. 
During td, irregularities larger than the typical noise level 
occur in m a /m in most of the experiments, see Fig. [SJ 
curve B. These irregularities are referred to as "micro- 
slips" in the following. Technically, a micro-slip, dm^, is 
defined as the increase of m a /m in time intervals where 
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FIG. 4: (a) The evolution of the normalized apparent mass 
m a /m as a function of time, (b) Closeup of one period. 
The total period At is the time between two peaks. t$ 
is the time through which the apparent mass decreases in 
one period, while t% is the time of increasing apparent mass. 
Aa — Am a /m is the amplitude of an oscillation. A subscript 
n is added to these parameters when they describe the specific 
values in oscillation number n. 

the time derivative of m a is positive. 

The observed oscillations in the apparent mass mea- 
sured under the ductile granular ensemble was seen for all 
cylinder diameters and filling heights, and proved very ro- 
bust to most perturbations applied to the system. Vary- 
ing the cylinder diameter and the filling height of grains 
did not affect the amplitudes and periods in any consis- 
tent manner. Amplitudes spanned 3%-24% of the total 
mass and the periods were in the range At = (0.7-47) 
h when all experiments are considered. Two otherwise 
equal experiments could produce different characteristics, 
in one case producing amplitudes of 6% and 20%, and pe- 
riods of 3.8 and 7.3 hours, respectively. The variability 
is probably due to details of the random packings that is 
beyond experimental control. 

Changing the surface properties on the cylinder wall 
was the only perturbation that dramatically affected the 
oscillations. Figure El shows results from experiments in 
which the surface friction was reduced by Teflon (curve 
A), and enhanced by (400 grit) sandpaper (curve C). In 
the following these experiments are referred to as 'the 
Teflon-' and 'the sandpaper experiment', respectively. No 
alteration was done to the surface of the wall in the ex- 
periment that produced curve B, which, apart from the 
surface, was identical to the Teflon- and sandpaper ex- 
periments. As can be seen from the figure, reducing or 
enhancing the wall friction both removed the weight oscil- 
lations. By reducing the friction on the wall the apparent 
mass increased slightly from the initial value (curve A, 
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FIG. 5: The resulting apparent masses for different surface 
treatments: Curve A was the result of coating the walls with 
Teflon (low friction). No coating of the Plexiglas wall resulted 
in curve B. Gluing sandpaper to the wall to enhance surface 
friction gave curve C. 

Fig. El- Although Teflon reduced friction considerably, 
it did not remove it fully, which would have made the 
apparent mass equal to the total mass of the grains. 

By increasing wall friction another behavior emerged, 
as the apparent mass decreased, apparently toward a con- 
stant level (curve C, Fig. [5J. Curve C was fitted excel- 
lently by 

m a /m = (m aoo + Am a exp [-t/r s ])/m , (1) 

where m aoo = (7.027 ± 0.001) g, Am a = (10.830 ± 0.005) 

g, and t s is a characteristic time constant of (13.52±0.01) 

h. The uncertainties are the likely error in the best fit 
parameters. Figure El shows the deviations between the 
data and the fit, (m a — m aoo + Am a exp[—t/T s ])/m. The 
exponential decay fits the observation exceptionally well, 
and the deviations are within the range [—0.0077, 0.0076] 
of the normalized data. Nevertheless, micro-slips are eas- 
ily recognizable above the experimental noise, which is of 
the order of 2 • 10~ 4 (0.03 g/142.6 g) of the normalized 
apparent mass, while the slips are of the order of 7- 10~ 3 . 
The experimental noise is not visible in the figure. 

A fit has also been made to the decreasing part td of 
each oscillation n in the curve of Fig. Ufa), which is 
consistent with logarithmic decay with time: 

^ = (l-B n ln[l + t/r d „]). (2) 

Here, m an {t) is the apparent mass of the n-th oscilla- 
tion, and m an (0), B n and Td n are best fit parameters 
to the equation, calculated for each oscillation n sepa- 
rately. TOan(0) is the best fit value of m an at the start 
of the decr ease, ba sed on the first 2.5 h of the decreas- 
ing m an . m an (0) — 76.4 [—0.4, 0.5] g is the median 
value of m an (0), with the quartile deviations in brackets. 
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FIG. 6: Deviations of the fit from the measured normalized 
apparent mass, ma/m, as a function of time t for the sand- 
paper experiment, see Eq. 
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FIG. 8: The deviations from the measured m a /m of its fit for 
the decreasing part of one oscillation, as a function of time, 
see Eq. H 
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FIG. 7: The decreasing apparent mass of 17 oscillations in 
one experiment (see Fig. Bla)) plotted as a function of time 
according to Eq. [5] (The expression on the horizontal axis is 
the time dependent part of Eq. 0) 



B n = 0.042 [-0.002, 0.004] is the median of the set of di- 
mensionless constants B ni and t^T = 0.16 [—0.02, 0.03] h 
is the median and quartiles of the set of Td„. Figured 
shows the collapse of the weight data when plotted ac- 
cording to Eq. [21 The limited dynamic range on both 
axes suggests that one can also fit the data by a power 
law with a small exponent. We have not found any theo- 
retical arguments for the choice of one fit over the other, 
thus the main observation is that the decreasing parts 
of the oscillations have the same form over the first 2.5 
hours, with a time constant of t\i~ — 0.16 h. The sandpa- 
per gave a decreasing exponential function with time, as 
gives the initial decrease of m a during t d in an oscillation: 



lim 1-Bln (1 + t/r d ) ~ l-Bt/r d ~ exp (-Bt/r d ) = exp (-t/r . 

The functional dependence is thus similar to the sandpa- 
per at the start of the decrease, with a time constant of 
t = r d /B = 3.8 h. 

The deviation from the fit is plotted for one oscillation 
in Fig. |H1 Large deviations on the order of 2% of the total 
mass (the micro-slips) develop some time into t d (typi- 
cally 3 hours in this experiment). All visible irregularities 
in this plot is above the noise level of the measurements. 

Taking the time derivative of m a as dm a /(mdt) = 
(m a {i + 1) — m a (i)) / m[t(i + 1) — t(i)], the set of positive 
increments of Tfi a jvri (the micro-slips, dui^^j and negative 
increments (dm~) were found for each oscillation's t d . 
The micro-slips were removed from the decreasing part 
of the oscillations by cumulative summation of dm~ , and 
the resulting data set fitted by a power law, 



£dm-(t)-l = -(t/T_) Q 



(4) 



The median and quartiles of the fitting parameters are 
a = 0.55 [-0.02, 0.02], and t_ = 130 [-37, 24] h, see Fig. 
No characteristic time exists for a power law, since 
-\ a (t/\T) a fits equally well for all A. 

The micro-slips dm^ were found as a function of time 
in all the oscillations of the experiment shown in Fig. 
and binned in 50 time intervals. The sum of micro-slips 
was taken for each bin, and divided by the size of the 
bin to produce the temporal evolution of micro-slip 'ac- 
tivity'. Figure presents the result. As the t d were 
of different lengths for each period, not all bins contain 
contributions from all oscillations. The bullets present 
times that include data from all oscillations, whereas a 
circle includes only data from t d long enough to con- 
tribute to the specific bin. The line through the data 
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FIG. 9: The cumulative sum of decreasing m a during an os- 
cillation as a function of the scaled time (t/r_) a , see Eq. [I] 
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FIG. 10: The micro-slip 'activity' as a function of time after 
each big slip. The activity is found as the sum of micro-slips, 
dmt /m, from all oscillations in Fig. [I] binned in times t n 
and normalized by the width of the bin. 



is a linear fit, based on all but the first bullet, given by 
E 9 "»+(tn)/t„ = A{t - t ). Here, A = (0.076 ± 0.005) 
h -2 , and to = (0.6 ± 0.2) h. The activity presented by 
the bullet at t ~ is probably remnants from the big 
slip that occurred at t = 0, thus the micro-slip activity is 
initiated at time to after each big slip and grows linearly 
until another big slip occurs. 

We could not find a model with few parameters that 
would fit the 'Teflon' results (curve A in Fig. |5J) due to 
the complex initial evolution. This curve also shows some 
micro-slips, of size 1.5 • 10 -3 , larger than the noise level 
of4-10~ 4 . 

The measurements of the height of the system as func- 
tion of time revealed that the vertical motion occurs in 
steps. This was seen in both strain experiments, and 
is shown in Fig. Illf a - ) and (c). Figures ITTT b) and (d) 
show the simultaneous measurements of the normalized 
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FIG. 11: Details of the strain as a function of time, mea- 
sured in two experiments compared to the weight m a . (a) 
The global strain, e, measured with 3 /im resolution (see Fig. 
|5J) as a function of time, (b) The normalized apparent mass, 
m a /m, as a function of time for the experiment in (a). (c)e 
measured with 46 resolution by the high resolution cam- 
era, (d) The apparent mass as a function of time for the 
experiment in (c). 



apparent mass, m a (t)/m. From the experiment with 3 
/im resolution, the minimum and maximum compaction 
velocities of the central part of the cylinder were found 
to be 5.4 • 10 -9 m/s and 7 • 10 -8 m/s, respectively. The 
maximum acceleration, occurring at the start of a com- 
paction step, was 1 • 10~ n m/s 2 . Comparing the region 
of decreasing m a of Fig. IllT b) to the strain in (a), a 
small but visible vertical movement occurs along the cen- 
tral axis of the packing during the weight decrease. The 
main increase of strain during one oscillation (that is, the 
step) takes place within the region in which the appar- 
ent mass increases from its minimum to its maximum. 
Unfortunately, the limited resolution of the strain mea- 
surements in Fig. IllT c) prevented a detailed comparison 
between the strain evolution of the 6 uppermost grains 
and the apparent mass. It is evident from this measure- 
ment, however, that the global strain motion is directly 
correlated with the changes in the apparent mass. A com- 
paction velocity of the uppermost grains of (0.6-3) TO -9 
m/s was found during t^, and 4-10 -8 m/s during ti. 

The dynamics of the wall contacts were studied in one 
experiment as described in section [H] Having found the 
'unwrapped', properly scaled surface of the cylinder (see 
Fig. |3J), we obtained a high contrast of the contact. The 
development of the area and center of mass position of 
each contact was followed through the experiment. 
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FIG. 12: The cumulative distribution P(s > Ay) as a function 
of normalized vertical contact displacement Ay/d c between 
two consecutive images. Motion downward has positive values 
of Ay. Curves A result during slips, while curves B present 
the remaining movement between time steps. The gray region 
through B covers the average value taken at each 1/79 interval 
of P, plus and minus the standard deviation from the average. 

The weight oscillations correlated strongly to the con- 
tacts' center of mass motion, while no such correlation 
was found with the changes in contact area. 

The contacts were seen to move "simultaneously" , that 
is, within the temporal resolution of the images, which 
means they all slipped within a period of 15-20 minutes, 
during tj. Figure IT21 shows the cumulative distribution 
P(s > Ay) of vertical contact displacement Ay between 
two consecutive images. The contact displacement is nor- 
malized by the average contact diameter d c . Each curve 
corresponds to the distribution in one time step of the 
experiment. Curves A present the motion during slips, 
while curves B are the motion in time steps between slips. 
The gray band through B spans the average ± the stan- 
dard deviation of vertical motion. 

The median vertical displacement of a contact during 
a slip was 6 [—1,2]% of the average contact diameter, 
d c . Outside of the slips the median displacement was 
only 0.07 [-0.20,0.24]% of d c . Figures Eta), (b) and (c) 
show the difference in one contact area between consec- 
utive images in one experiment. White corresponds to 
new contact area, black to area that was left since the 
previous image, and light gray shows contact area where 
no changes occurred. Figures I13f d) and (e) show the 
normalized apparent mass and the average strain of the 
upper (diamonds) and lower (circles) 20 wall contacts, re- 
spectively. The markers in both plots represent the times 
when pictures were taken. In Fig. Il^f dl circles mark the 
times when 15% of the contacts moved more than 1% 
of the average contact diameter in 20 minutes (since the 
last image) . The bullets show the times when 80% of the 
contacts moved at least 2% of the average contact diam- 
eter. Triangles represent the times when pictures were 




30 35 40 45 
t[h] 

FIG. 13: (a), (b) and (c) shows difference images of a contact 
between consecutive images. White is newly established area, 
black is area that no longer is a part of the contact, and 
light gray is the unchanged contact area. The center of mass 
motion between the images are shown as black lines. The 
length of the lines is exaggerated 10 times, (d) shows the 
normalized apparent mass, the triangles (A) mark the times 
when pictures were taken of the ensemble. Circles (o) mark 
the times when minimum 15% of the contacts moved more 
than 1% of the average contact diameter. Bullets (•) mark 
the times when more than 80% of the contacts moved at least 
2% of the average contact diameter. The lower plot (e) shows 
the average strain development found from image analysis for 
the 20 lower (o, curve B) and upper (o, curve A) wall contacts. 
Filled symbols represent the times that a picture was taken at 
or immediately after a peak in the apparent mass presented 
in (d). 



taken. Based on the observed area of the grain-wall con- 
tacts and the measured m a , the average load per square 
millimeter carried by a contact was calculated to be in 
the range (0.5-1.2) kPa. 

Table 13 presents the characteristic values of various pa- 
rameters: (a) gives the median period, amplitude, inter- 
vals td and ti, and characteristic times Td and t for the 
oscillations in one experiment (see Fig. 0J. (b) presents 
the characteristic time from the fit of m a /m from the 
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Characteristic values 





At 


6.4 [-0.7, 1.2] h 




td 


5.2 [-0.3, 1.1] h 




U 


0.8 [-0.3, 1.0] h 


(a) 




0.16 [-0.02, 0.03] h 




TO = Td/ B 


~ 3.8 h 




to 


0.6 ± 0.2 h 




Aa 


12.6 [-0.3, 1.0] % 


(b) 


T a 
T e 


13.52 ± 0.01 h 
~ 10" 6 h 




Is 


~ 260 


(c) 


/(, 


~ 74 ^m 


h 


~ 101 (im 




k 


~ 115-200 



TABLE I: (a) Median values of the period At, amplitude 
Aa, the intervals td and ti, the characteristic times Td and 
to = Td/ B of decreasing m a , and to of activation of micro-slips 
of the experiment presented in Fig. Ufa), (b) Characteristic 
times t 3 of the m a /m evolution in the sandpaper experiment, 
and r e , the estimated time of relaxation of elastic stress, (c) 
Estimated characteristic length scales, from time scales in (a) 
and (b), see section ITV1 



sandpaper experiment, and the estimated characteristic 
time of elastic relaxation (see section HV|) . 

One experiment was performed to understand how the 
granular geometry of the ensemble affected the apparent 
mass. The granular ensemble was exchanged with a non- 
porous slab of Play-Doh that did not fill the cylinder, 
but touched both the bottom and the walls of the setup. 
This experiment is referred to as 'the bulk experiment' in 
the following. The slab was left to flow into the available 
space, and the apparent mass was measured as before, 
see curve B of Fig. ^1 Again, a granular version of this 
experiment was conducted for comparison, in which the 
total mass of the grains and the cylinder diameter were 
the same as those of the bulk experiment, see Fig. 1141 
curve A. As seen from the figure, both setups produced 
weight oscillations, thus the granular geometry is not the 
(only) reason for the oscillations. Oscillations started 
later in the bulk case than in the granular case, and both 
systems show uncommon irregularities in their periods. 
The granular system had nearly 100 grain-wall contacts, 
while the bulk sample had 3-4 large contact areas. The 
oscillations are probably due to the multi-contact nature 
of the interface between the deforming sample and the 
confining cylindrical wall. 



IV. DISCUSSION 

The self-compaction of a ductile ensemble depends on 
the deformability of the grains and on a porous struc- 
ture. The granular geometry of the ensemble was not 
necessary for oscillations to form, as weight oscillations 
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FIG. 14: The normalized apparent mass, m a /m, as a function 
of time, t, at the bottom of a cylindrical Play-Doh sample 
(B), as compared to m a /m from a granular geometry (A) as 
a function of time. 



also resulted under a bulk slab of material that deformed 
viscously into available space. This result emphasizes 
the importance of the multi-contact wall interaction to 
the observed oscillations in the apparent mass. 

The grain-wall interaction proved to be crucial to the 
oscillations in the apparent mass by the experiments with 
varying wall friction. No oscillations were observed when 
increasing or decreasing the wall friction from that of the 
regular experiments with Plexiglas walls. The evolution 
of m a in these experiments is interesting because it shows 
two different behaviors according to the wall friction. A 
low wall friction resulted in an increasing apparent mass, 
while a high wall friction made the measured weight de- 
crease. The same mechanisms leading to these results 
are likely to be the reason for the oscillations observed in 
m a in the regular experiments. 

The reason for the decrease of the apparent mass must 
be that the walls sustain an increasing part of the total 
grain mass, that is, a dynamic vertical force must act 
on the grain contacts from the wall. This force could be 
friction, shear cohesion, or a combination of the two, and 
will be referred to as the grain-wall interaction in the 
following. 

The increasing weight was initially believed to be due 
to a new internal grain-grain contact. As the stress dis- 
tribution in hard granular media is known to be very sen- 
sitive to local arrangements of grains, a new contact was 
believed to change the stress distribution. New contacts 
would preferentially form in the vertical direction, be- 
cause of the anisotropic compaction, and thus would tend 
to redirect stresses toward the bottom. The number of 
new contacts in the ensemble is limited, and the average 
number of contacts per grain increased from 6.5 to 7 [20| 
during a typical experiment. 50-100 new contacts were 
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expected to form during an experiment, which is roughly 
twice the typical number of oscillations. If we assume 
that not all new contact formations are noticed in m a , 
perhaps because of their position in a low stressed region, 
that would explain the shortage of oscillations, and the 
micro-slips in m a in the oscillations and the sandpaper 
experiment. On the other hand, this assumption directly 
disagrees with the nearly constant amplitudes seen in all 
experiments. The experiment of a bulk slab, which also 
produced weight oscillations, eventually proved that new 
internal contacts between grains were not the main rea- 
son for the weight oscillations. Stress redistribution is, 
however, thought to take place continuously during the 
slow internal flow of material, both in the granular and 
the bulk systems. 

In principle, elastic energy could be stored in com- 
pressed parts of the packing after a slip, resulting in a 
decreased grain-wall interaction. The relaxation of this 
elastic energy could cause the observed decrease in the 
apparent mass. The characteristic time of elastic relax- 
ation is expressed as the ratio of viscosity fi to bulk mod- 
ulus K . We know that the viscosity of Play-Doh is of the 
order of 10 5 Pas for shear rates as low as 10 ~ 6 s _1 . The 
bulk modulus was not measured, but it is expected to be 
closer to that typical of fluids (Kj ~ 1-2 GPa, j^) than 
that of iron (iQ ~ 170 GPa, (23), thus on the order of 
10 9 Pa. The resulting estimate of elastic relaxation time 
for Play-Doh is 

Te = fj,/K = 10 5 /10 9 = l(T 4 s . (5) 

Elastic compressive stresses should relax in (less than) 
seconds, which is much less than any time scale observed 
in an oscillation. 

Another explanation for the decreasing m a emerges 
from the assumption that the ratio of horizontal to verti- 
cal stresses increases with increasing packing fraction. If 
friction is assumed to be proportional to the normal force, 
an increasing horizontal stress in the packing would result 
in increased wall friction, hence a decrease in m a . The 
packing fraction increases approximately 10% during the 
experiment, while the characteristics of the oscillations 
does not change. This implies that the packing fraction 
is not the main parameter for describing the dynamic 
behavior. 

The reason for a decreasing apparent mass can be seen 
in connection to the shearing of grain-wall contact re- 
gions. During the time t& of decreasing m a the strain in- 
creases very slowly, suggesting that only an internal flow 
of grains contributes to the strain in this regime (see Fig. 
Illf a) and (b)). The analysis of the motion of grain-wall 
contacts shows that the vertical motion of contacts in 
this regime is limited and noisy, thus most contacts are 
practically at rest (even though the central part of the 
packing creeps). Due to the slow flow internally in the 
packing, they are also continuously sheared. 

There are clear slips of the order of 2% in the nor- 
malized apparent mass during the decreasing part of the 
period in the granular setups (see Figs. 0]and|SJ|. Micro- 



slips are not seen in the weight data from the bulk exper- 
iment, thus their origin seems to be the granular geome- 
try, or possibly the large difference in the number of wall 
contacts between the bulk and granular systems. No col- 
lective motion is seen at the wall during the micro-slips 
in the granular experiment, although 5% of the contacts 
move a distance of 1% of d c in every time step, thus their 
motion might be connected to the measured micro-slips 
in m a . Micro-slips might be due to the internal reorga- 
nization of forces within the granular system, which may 
trigger some of the grain-wall contacts into motion. A re- 
organization of forces must also take place in the material 
in the bulk experiment, although probably in a different 
way than that of the more complex granular geometry. 
The reorganization must increase the average shear stress 
in the contact regions, which again leads to an increase 
of the vertical grain-wall interaction. Once a contact ex- 
periences a shear stress that can not be sustained by the 
grain-wall interaction, it "breaks", or starts to move. 

The strain development could not be measured in the 
sandpaper experiment, thus whether this system com- 
pacted much is not known. Similar, but smaller, micro- 
slips than those seen in regular experiments were seen in 
the sandpaper experiment. This suggests that internal 
stress rearrangement was taking place. The grain-wall 
interaction was considerably higher in the sandpaper ex- 
periment than in the regular setup (as the apparent mass 
reached a minimum of 15% of the total mass). It is rea- 
sonable to assume that the contacts did not move much, 
or in any correlated manner, based on the lack of weight 
oscillations. 

The direct correspondence between the step in the 
strain and the increasing m a in Fig. Illf a) and (b) im- 
plies that the motion of wall contacts is very important 
for the weight increase. Assuming that wall contacts are 
broken, or mobilized, at a critical shear stress, one or 
more contacts will initiate a slip, and the others follow. 
The contacts that break contribute to a decrease in the 
total wall interaction, thus a possible increase of the ap- 
parent mass. The sum of wall interactions decrease over 
a time period that must depend on how fast the contact 
breaking propagates among the other contacts, and how 
fast each contact breaks. From our temporal resolution 
in the study of grain-wall contacts, we see that all grains 
move within 20 minutes in connection with the slip. The 
apparent mass will increase according to the decreasing 
wall interaction F w , as the force balance of the system is 
Y,F = m a g~ F w = ma < 10~ 12 N ~ 0. 

The strain development was not measured in the Teflon 
experiment, thus it is not known whether the strain had 
similar steps during the compaction as in the regular 
experiments. Based on the direct correlation between 
weight oscillations and the observed strain in the regu- 
lar experiments, however, it seems likely that the wall 
contacts in the Teflon experiment in some sense moved 
continuously, as no oscillations in m a were observed here. 
Micro-slips were observed in m a , however, thus some dy- 
namic interaction between the grains and the wall was 
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present, probably because of internal stress rearrange- 
ments. Sliding contacts also support some of the grain 
mass, as neither during tj in the regular experiments nor 
in the Teflon experiment does the apparent mass reach 
100% of the grain mass. The grain-wall interactions dur- 
ing motion are smaller, however, than in the static case, 
as the apparent mass increases during motion in the reg- 
ular experiments, see Fig. 1111 

That all contacts are mobilized within a time inter- 
val corresponding to a slip could imply that, when suf- 
ficiently sheared, they are sensitive to changes in the 
stress distribution, and thus easily triggered. From Fig. 
I13f d) we see that more than 80% of the contacts move 
more than 2% of d c during a slip event, and that 15% 
move at least 1% of d c immediately before these slips. In 
some cases, although not consistently, 15% of the con- 
tacts move at least 1% of d c in connection to micro-slips. 
Also, the activity of micro-slips increases during ta, which 
suggests that the system becomes more critical. 

The time scales of the system spans a factor 100, see 
Table[U ranging from 0.16 h to 13.52 h. It is tempting to 
speculate that these time scales reflect the spatial dimen- 
sions in the system, from 1 mm (diameter of small contact 
area) to 10 cm (filling height). A direct estimate of the 
maximum length scale can be made from the velocities 
and the observed time scales. Assuming that the grain- 
wall contacts in the sandpaper experiment do not slip, 
the internal flow of velocity Vd = 5.4 • 10~ 9 m/s with the 
characteristic time t s gives a length scale l s — 260 //m. 
The corresponding length from the initial exponential de- 
crease of m a in an oscillation is Iq — tq- Vd m/s = 74 /mi, 
and from the t& and tj, we get Id — td ■ Vd — 101 /im and 
h — U ■ Vi =115-200 /im, respectively. Uj is the velocity 
of the bulk during a slip. The range of k results from the 
different compaction velocities found during ti in the two 
experiments presented in Fig. 1111 The length scales ex- 
tracted from the characteristic times span a smaller range 
than the time scales do, and are much smaller than the 
macroscopic lengths mentioned above. The small length 
scales suggest that details of the contact motion might be 
of importance to the time scales observed in the system. 

Flow of viscous fluid along a wall can be described by 
a Navier length [23. An average contact velocity, v c , 
during a slip can be found from knowing that contacts 
slip 6% of the average contact diameter in 20 minutes, 
v c = 1.15 • 10~ 7 m/s. The amount of fluid slip along a 
wall is given by the Navier length, b — n/k, where mu is 
the fluid viscosity and k is the surface friction coefficient 
given by a/v c - The average shear, er, of a contact was 
found to be between 0.5-1.2 kPa, thus k is within the 
range (2.7-ll)T0 9 Pas/m. The Navier length is then 
b € [27-90]/im, slightly smaller, but of the same order as 
some of the lengths estimated above. 

The motion of a contact was not studied with sufficient 
temporal or spatial resolution to conclude whether the 
whole contact slid a fraction of d r at constant velocity, or 
it slid by self-healing slip pulses [23, |24j . Both processes 
are known from experiments on frictional motion in low 



velocity sheared systems of (hard) granular systems pj| 
and slipping of a gel/glass interface |24| . 



V. CONCLUSIONS 

We observe semi-regular oscillations in the measured 
apparent mass, m a , at the bottom of a self-compacting 
ductile grain packing. The oscillations in one particular 
experiment are on the order of 10% of the total mass 
to of the grains, and have periods of roughly 6 hours. 
The oscillations persist when the granular setup is ex- 
changed with a bulk sample of the same ductile material, 
but disappear when the grain-wall interaction is reduced 
or increased. Grain-wall contacts are seen to move col- 
lectively in correspondence to the slip events in m a , as at 
least 80% of the contacts move a distance larger than 2% 
of the average contact diameter during a slip, see Figs. 
dandED 

The decrease of the apparent mass in an oscillation is 
thought to be the result of shearing of static wall contacts 
between grains and the container wall. The slow ductile 
flow internally in the cylinder causes a dynamic stress 
distribution, which results in a continuous increase of the 
shear stress at the grain-wall contacts. This continuous 
increase is the reason for the decreasing apparent mass. 

"Micro-slips" of the order of 2% are seen in the nor- 
malized apparent mass m a /m during the decrease, which 
probably result from internal stress redistribution in 
granular setups, as they were not seen in m a of the bulk 
experiment. The micro-slips correspond in some cases to 
limited grain- wall contact motion, and their 'activity' in- 
creases during the interval of decreasing m a . These slips 
are also seen when the grain-wall interaction is reduced 
or enhanced, that is, when contact motion is stimulated 
or repressed. 

Different characteristic times have been found from 
curve fitting of the apparent mass evolution during the 
'sandpaper' experiment and the decreasing part of os- 
cillations in one experiment. We have also estimated 
a typical timescale of relaxation of elastic compressive 
stresses, and concluded that elasticity is not the driv- 
ing mechanism for the observed oscillations. The charac- 
teristic times, together with the period and intervals of 
increasing and decreasing m a , are presented in Table [I] 
A successful model should reproduce these characteristic 
times. 

Some attempts at constructing a minimum model have 
been pursued, but the models were discarded as they de- 
pended on a finite acceleration or on unknown micro- 
scopic parameters of the system. Further work is neces- 
sary to understand the dynamic behavior of the system, 
and ideas and modeling approaches are welcome. 
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